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ABSTRACT 

Context. B[e] supergiants are surrounded by large amounts of hydrogen neutral material, traced by the emission in the optical 
[Oi] lines. This neutral material is most plausibly located within their dense, cool circumstellar disks, which are formed from the 
(probably non-spherically symmetric) wind material released by the star. Neither the formation mechanism nor the resulting structure 
and internal kinematics of these disks (or disk-like outflows) are well known. However, rapid rotation, lifting the material from the 
equatorial surface region, seems to play a fundamental role. 

Aims. The B[e] supergiant LHA 1 15-S 65 (in short: S 65) in the Small Magellanic Cloud is one of the two most rapidly rotating B[e] 
stars known. Its almost edge-on orientation allows a detailed kinematical study of its optically thin forbidden emission lines. With a 
focus on the rather strong [Oi] lines, we intend to test the two plausible disk scenarios: the outflowing and the Keplerian rotating disk. 
Methods. Based on high- and low-resolution optical spectra, we investigate the density and temperature structure in those disk regions 
that are traced by the [Oi] emission to constrain the disk sizes and mass fluxes needed to explain the observed [Oi] line luminosities. 
In addition, we compute the emerging line profiles expected for either an outflowing disk or a Keplerian rotating disk, which can 
directly be compared to the observed profiles. 

Results. Both disk scenarios deliver reasonably good fits to the line luminosities and profiles of the [Oi] lines. Nevertheless, the 
Keplerian disk model seems to be the more realistic one, because it also agrees with the kinematics derived from the large number of 
additional lines in the spectrum. As additional support for the presence of a high-density, gaseous disk, the spectrum shows two very 
intense and clearly double-peaked [Can] lines. We discuss a possible disk-formation mechanism, and similarities between S 65 and 
the group of Luminous Blue Variables. 
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• 1. Introduction drogen right from (or at least very close to) the stellar surface 

•n (Kraus & Borges Fernandes 2005 ; Kraus et al. l2006l l2007l l. so 

O . B[e] supergiants in the Large and Small Magellanic Clouds jhat molecules and dust are forming in the vicinity of the star. 

O . (LMC and SMC), even though studied in great detail, ai-e still indeed, for the LMC B[e] supergiant R 126, Kastner et al. (200i 

7^, ■ far from being understood. A non-spherically symmetric wind is fou^d that the inner edge of the massive dusty disk must be lo- 

> , suggested by, e.g., polarimetric observations (Magalhaes 1992; ^ated at ~ 360R„ which is about three times closer to the star 

Magalhaes et al. 2006; Melgarejo et al. 2001), and a circum- jj^an the value of ~ 1000 7?. formerly suggested by Zickgraf et 

_ stellar disk seems to be confirmed by their strong infrared ex- ^1. ( 1985). The much closer inner edge of the dust disk requests 

^ cess e mission due to circumstellar dust (e.g., Zickgraf et al. that the [Oi] emission must originate from distances between the 

■ - . [1986) as well as by the emission of intense molecular bands like surfac e of R 126 and ~ 360 R„ which was confirmed by Kraus 

the CO first-overtone bands in the near infrared K-h and sp ectra gt al. (|2007] | considering a hydrogen neutral disk right from the 

(McGregor et al. 1988a 1988b, 1989; Morris et al.fT996l l. and stellar surface. 

TiO bands at optical wavelengths (Zickgraf et al. |1989] l. Based _ , , ,. , • , , , • 

^. 1 • f f- 1 rr^nr » u ti u To guarantee that the disk material close to a luminous and 

on the analysis of optical [Oil lines, it has recently been sug- , ^ . . ,,,-,, , • / 

. ■ ., . ;^, ■• , J TDr 1 • » . 1 • u hot supergiant star is neutral, the disk has to be massive (e.g., 

gested that the disks around Bie] supergiants are neutral in hy- ,, ^„ ° ,^„^^. „ ' , ,. , , . ° 

Kraus & Lamers 120031 1. Because the disks around supergiants 

cannot be pre-main sequence in origin, their formation must 
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* Based on observations done with the 1.52m and 2.2m telescope at be linked to non-spherical (predominantly equatorial), high- 

the European Southern Observatory (La Silla, Chile) under the agree- ^jg^^j^y ^^^^ j^^^ ^^.^^ ^j^g ^g^^^^j ^^^.^ gj^j^gj. ^y some contin- 

ment with the Observato^rio Nacional-MCT (Brazil) and under pro- ^^^^ ^^^^^ ^^^^^-^^ ^^^^ ^^ ^ ^^^^ mass-ejection event(s). 

grams 076.D-0609(A) and 080.A-9200(A) -' -' ■' ^ 

** We have to report with great sadness that our friend and colleague. The formation of this anisotropic, high-density mass loss is 

Francisco Xavier de Araiijo, deceased before this work could be fin- assumed to be linke d to r apid stellar rotation (e.g. Maeder & 

ished Mevnet i2000i Kraus |2006] |. Especially rotation close to the crit- 
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Table 1. Parameters of S 65. The rotation velocity was derived 
by Zickgraf (120001) . the remaining parameters by Zickgraf et al. 
(IT986] l. 



T-eff 

[K] 


[Re] 


log L./Lo 


osin; 
[kms-'] 


E{B - V) 


inclination 


17 000 


81 


5.7 


~ 150 


0.18 ±0.02 


± edge-on 



width of 4". The instrumental setup employed provides a res- 
olution of ~ 4.6 A in the range of 3800-8700 A. The S/N ratio, 
in the 5500 A continuum region, is aproximately 60. The B&C 
spectrum was reduced with standard IRAFQ tasks, such as bias 
subtraction, flat-field normalization, and wavelength calibration. 
The absolute flux calibration was made using spectrophotomet- 
ric standard stars cited in Hamuy et al. (|1994| l. 



ical velocity could trigger the (rotationally induced) bi-stability 
mechanism (Lamers & Pauldrach 119911 Pelupessy et al. 120001 
Cure 120041 Cure et al. 2005 ) as a plausible disk-formation sce- 
nario. And indeed, for at least four members of the B[e] super- 
giant group, the projected rotation speed (i.e., ysin/) could be 
determined, showing that two of them are rotating at a substan- 
tial fraction (i.e., at least 75%) of their critical velocity. These 
are the SMC stars LHA 1 15-S 23 discussed in detail by Kraus et 
al. (2008) and LHA 1 15-S 65 (Zickgraf 2000). Whether all B[e] 
supergiants are rapidly rotating, is unknown however. The high- 
density wind observed around the other B[e] supergiants hides 
the central star, and the circumstellar material additionally con- 
tributes and pollutes the optical spectrum with a huge amount of 
emission lines. Consequently, the detection of uncontaminated 
photospheric lines, which would be appropriate for the determi- 
nation of the projected stellar rotation velocity, is strongly ham- 
pered. 

We study here the neutral material around the rapidly rotat- 
ing SMC B[e] star LHA 1 15-S 65 (in short: S 65). The star S 65 is 
one of the four confirmed B[e] supergiants in the SMC. Its stellar 
parameters are summarized in Table[T]and have been determined 
by Zickgraf et al. (1986) and Zickgraf (2000). Our investiga- 
tion makes use of the three optical [Oi] emission lines detected 
in our spectra. We particularly aim to simultaneously reproduce 
the line luminosities and line profiles, to test the validity of any 
of the two most plausible B[e] supergiant disk scenarios: (i) an 
outflowing disk, which is neutral in hydrogen straight from the 
stellar surface, and (ii) a Keplerian rotating disk. 



2. Observations and reductions 

We obtained high-resolution optical spectra using the Fiber-fed 
Extended Range Optical Spectrograph (FEROS), and also low- 
resolution optical spectra using the Boiler & Chivens spectro- 
graph (B&C). 

The FEROS spectra were obtained in three different years. 
The oldest one is from 1999 October 27, when the spectro- 
graph was still attached to the 1.52-m telescope at the European 
Southern Observatory in La Silla (Chile), while the other spec- 
tra were taken on 2005 December 12, and 2007 October 3 and 
4, after FEROS was moved to the 2.2-m telescope. FEROS is 
a bench-mounted Echelle spectrograph with fibers, where each 
one covers a sky area of 2" of diameter The wavelength cover- 
age goes from 3600 A to 9200 A and the spectral resolution is 
R - 55 000 (in the region around 6000 A). The spectra were ob- 
tained with exposure times of 7200 (1999), 900 (2005), and 900 
(2007) seconds respectively. In 2005 and 2007, two exposures 
have been taken, which were then added for a better S/N ratio. 
The S/N ratio in the 5500 A region is approximately 90 (1999), 
40 (2005), and 65 (2007). We adopted the complete automatic 
on-line reduction. 

The low-resolution B&C spectrum was taken on 1999 
October 28, with an exposure time of 900 seconds and a slit 



3. Results 

3.1. Spectral characteristics 

A first detailed list of identified lines in the optical spectral 
range of S 65 has been compiled by Muratorio ( 1981 ). Later on, 
Zickgraf et al. ( 119861 ) re-observed S 65 at higher spectral resolu- 
tion and with higher signal-to-noise ratio, and described their 
optical spectra in great detail. They found that the line spec- 
trum showed many features in common with shell stars. These 
are the Balmer lines, which display a PCygni profile of Beals' 
type III with a shell absorption reaching in most cases below 
the continuum, the Feii emission lines with a central absorp- 
tion core, and many additional singly ionized metal lines (like 
those from Cm and Tin, see Fig.lTJ showing very narrow cen- 
tral absorption components which are nearly unshifted with re- 
spect to the emission component. For the typical B[e] super- 
giants, these well pronounced shell lines were thus interpreted 
by Zickgraf et al. ( 1986) as clear signatures for an edge-on orien- 
tation of the S 65 system. Further support for an edge-on orien- 
tation came from optical polarization measurements performed 
by Magalhaes (119921 ) and Melgarejo et al. (I2001I ). Follow-up ob- 
servations with FEROS (Stahl l200T]) showed no variability in the 
spectral lines of S 65, when compared to the spectra of earlier in- 
vestigations. Indeed, comparing our own three spectra that cover 
a period of eight years to those of Zickgraf et al. tl98 6), as well 
as to the FEROS spectrum of Stahl (120011 ) that was obtained one 
year prior to our first spectrum, we could not find any significant 
difference in the individual line profiles, indicating that S 65 is 
not undergoing any severe variations. Nevertheless, before going 
into any modeling details, we will first provide a short descrip- 
tion of the elements seen in our spectra and their line profiles. 
Due to the good quality of the spectra taken in 2007, we mea- 
sured all the lines reported in this paper from these spectra. The 
spectra taken in 1999 and 2005 are used for a qualitative com- 
parison. 

Figure[T] covers some small portions of the optical spectrum 
of S65. The numerous lines seen in these short wavelength 
ranges represent the richness in emission lines of the complete 
spectrum. A detailed investigation revealed only a few very weak 
photospheric absorption lines. Ten lines are from Hei, which are 
all at least partly contaminated by or blended with circumstellar 
emission lines. This is shown for eight of them in Fig.|2] We also 
tentatively detected two lines from Nn, AA 3995 and 5679, and 
possibly one line from On at A 3954. 

The Balmer lines are present from Ha up to H15. Figure[3] 
shows the line profiles of Ha to H6. The lowest Balmer lines 
display double-peaked profiles with a so-called violet-to-red ra- 
tio V/R < 1 and a broad central absorption that extends below 
the continuum (except for Ha, see Fig.[3]). These line profiles 
are well-known from classical Be stars (or Be shell stars, e.g.. 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. Parts of our FEROS spectrum showing the two types of emission lines discussed in the text, i.e. the single-peaked emission 
lines of the forbidden transitions, and the permitted emission lines which all display a central or slightly blue-shifted sharp absorption 
component. 



Hanuschik I2000I I viewed edge-on. From \\6 on the shape of 
the lines turn into a PCygni profile. The strength of the emis- 
sion component thereby quickly decreases with higher quantum 
numbers, until for Balmer lines higher than HI 3 only the absorp- 
tion component remains. The absorption thereby consists of the 
P Cygni part and a shallow, but rather broad photospheric part. 



as can be seen akeady for the W6 line in Fig. [3] Also present are 
lines of the higher Paschen series, i.e., PalO up to Pa24, display- 
ing P Cygni-type profiles. For the Paschen lines the absorption 
component also consists of a rather weak P Cygni absorption 
component, extending into the very shallow but much broader 
photospheric component. 
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Fig. 2. Photospheric Hei absorption lines. The lines are cen- 
tered on the meanjystemic velocity, which was determined as Fig.3. Balmer line profiles showing a variation in shape. The 

lines are centered on the mean systemic velocity of S 65, which 
was determined as described in Sect. l3.2l 



described in Sect. 13. 21 



The huge number of additional lines seen in our spectra can 
basically be divided into two groups: single-peaked emission 
lines, and emission lines with a strong but very narrow, almost 
central absorption component reaching in most cases far below 
the continuum. The two types of profiles are nicely seen in the 
parts of the spectrum shown in Fig.[T] 

All single-peaked emission lines are found to come from for- 
bidden transitions. Most of the lines belong to [Fen], but a few 
lines of [Cm] and [Niii], as well as one line from [Sn] are also 
visible. Additional lines (but blended or very weak) from [Tii], 
[Mgi], [Mgii], and [Nii] might be present. We also found a few 
very weak emission lines from permitted transitions; they belong 
to Ci and Ni. 

The group of emission lines with the strong central absorp- 
tion contains by far the largest number of lines. We could iden- 
tify more than 230 lines, most of them belonging to transitions 
of Pen, Cm, and Tin. But few lines from other elements like 
Mnii, Mgi, Can, Sin, and Fei were also found. The strong, sharp, 
and almost central absorption components indicate high-density 
absorbing material along the line of sight with a rather small ra- 
dial velocity component. As mentioned earlier, the line shapes 
recall those seen in Be shell stars or Be disks viewed edge-on. 
It is therefore reasonable to assume that S 65 is surrounded by a 
high-density disk or ring structure viewed edge-on, as suggested 
bvZickgrafetal. C 19861) . 

Interestingly, the spectrum shows prominent emission of the 
[Oi] lines as well as of the [Can] lines. However, their line pro- 
files do not fit into the groups described above. Instead, their 



profiles show a clear double-peaked structure (see Fig.|4|l. This 
type of profile is usually interpreted with Keplerian rotation in a 
circumstellar disk. 



3.2. Systemic and projected rotational velocity of S 65 

The strong pollution of the photospheric Hei absorption lines 
by adjacent emission from the circumstellar material (Fig.|2]) 
makes it impossible to accurately determine their line centers. 
Therefore we need to rely on the emission lines to derive the 
systemic velocity. As mentioned in the previous section, the per- 
mitted lines all show a sharp central absorption. The slight asym- 
metry seen in most of these line profiles in form of a slightly 
stronger red emission compared to the blue emission peak ham- 
pers a proper determination of their line centers. In addition, 
this asymmetry indicates that the central absorption components 
seem to be in many cases slightly blue-shifted. Hence they do 
not mirror the systemic velocity of the system, but only the dy- 
namics along the line-of-sight through the high-density disk (see 
Sect.|4]i. The only reasonable set of lines for a proper systemic 
velocity determination is thus provided by the optically thin, 
single-peaked and symmetric profiles of the forbidden emission 
lines. From the whole sample of forbidden lines detected, we re- 
stricted our selection to the strongest and unblended ones and 
measured their central wavelengths by fitting Gaussian profiles 
using IRAF From the difference between the measured and the 
laboratory wavelengths we computed the radial velocity, which 
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Fig. 4. Line profiles of the [Can] (top) and [Oi] lines (bottom). 
The lines are centered on the mean systemic velocity of S 65, 
which was determined as described in Sect. l3.2l 



is plotted in Fig.|5] The mean radial (heliocentric) velocity we 
found is V - (+ 189.23 +0.74) km s"', which can be interpreted as 
the systemic velocity of S 65. Our mean value agrees very well 
with the value of (+191 ± 4)kms"' derived by Zickgraf et al. 
( 119861 1. and also fairly well with the values of +184kms ', and 
(+193+6) km s"' obtained in earlier investigations by Feast et al. 
(1196 0) and Maurice (il976i) . respectively, when we account for 
the lower resolutions and lower S/N values of those older spec- 
tra. Therefore we consider our value as reasonably good and use 
it throughout the paper to correct the measured velocities with 
respect to the systemic velocity, as done, e.g., for the line profile 
plots in Figs.|2]to|4] 

We emphasize that the Hei absorption lines cannot be formed 
in the (polar) wind of S 65, because the effective temperature 
of the star (and hence the wind temperature) is too low to 
keep helium ionized beyond the photoshpere. Therefore they 
can be considered as of pure photospheric origin. The width of 
these Hei lines and also of the (shallow) photospheric absorption 
components of the Balmer and Paschen lines cannot be caused 
by broadening mechanisms like macroturbulence and pressure 
broadening alone. These two effects in the atmospheres of B- 
type supergiants typically add up to a few tens of kms"', only. 
Instead, the broad and roundish shape of the lines suggests that 
they are predominantely broadened by the possibly high rotation 
of the star, as suggested by Zickgraf (i20001 l. 
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Fig. 5. Radial (heliocentric) velocity of S 65 derived from the 
strong and symmetric non-blended single-peaked forbidden 
emission lines. 



Table 2. Measured HWHM and derived u sin / values of the three 
Hei lines in Fig.|6] 



Line 



HWHM 

[A] 



vsmi 
[kms-'] 



Hei 4922 L81+0.10 141.8 + 7.4 
Hei 5876 2.37 ± 0.09 155.2 ± 6.0 
Hei 7065 2.91+0.14 158.1+7.5 



Our derived systemic velocity allows us to determine the 
proper line centers of the Hei lines. Inspection of these lines 
(see Fig.|2]i shows that the blue parts of at least three lines (Hei 
AA 4922, 5876, and 7065) seem to be uncontaminated, nicely 
displaying the roundish profile shape as expected in rotation- 
ally broadened photospheric absorption lines. The blue half of 
these three Hei lines are shown in Fig.|6l Their wings extend 
to Doppler velocities of about 150-180kms"'. A comparable 
Doppler velocity is also seen in the shallow photospheric ab- 
sorption components of the hydrogen lines (see, e.g., the line 
profile of H(J in Fig.[3]l. 

We measured the half-width at half-maximum (HWHM) 
values of these three Hei lines in Fig.|6l From those values, the 
projected rotation velocity, u sin /, can be computed under the as- 
sumption of a stellar rotation profile as derived, e.g., by Gray 
(119761 ) and with a limb-darkening coefficient of 0.6. The mea- 
sured HWHM and the derived v sin / values are listed in Table|2l 
Compared to the value of ~ 150 km s' obtained by Zickgraf 
(I2000I I. our values excellently agree. Given that the Hei /15876 
line is the strongest in our sample, we consider our value of 
(155.2 + 6.0)kms"' derived from that line as a reasonably good 
value for the projected rotation velocity of S 65. 

We note however that this value might be un underestimation 
of the real projected stellar rotation velocity. Expressed in terms 
of the critical rotation velocity, the derived value of ~ 155 km s"' 
corresponds to vsini/ua-it - 0.75. For stars rotating so rapidly 
Townsend et al. (120041 1 have shown that the HWHM method is 
not reliable anymore to derive the correct projected stellar rota- 
tion velocity, because the effects of gravitational darkening in- 
fluence the photospheric absorption lines in a way that their line 
width does not increase anymore for usini/uait > 0.75. This 
means that even if the star would rotate at a velocity close to the 
critical one, the measured HWHM would only imply a value of 
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Table 3. Velocities obtained from the peak separation and 
FWHM values and extracted line luminosities of the [Oi] lines. 
The errors in the line luminosities are on the order of 10%. 



Fig. 6. Uncontaminated blue parts of three Hei lines from Fig.|2] 



vsmi/u„it - 0.75. Our value can thus only be considered as a 
lower limit to the real projected stellar rotation of S 65. 

3.3. The [Oi] lines 

In the remaining part of this section we concentrate on the [Oi] 
emission lines arising in the circumstellar material of S 65, aim- 
ing to find a resonable scenario to explain the observed [Oi] line 
luminosities and profiles. 

The FEROS spectrum displays three [Oi] lines (see Fig.|4|. 
They originate from transitions within the five lowest energy 
levels of neutral oxygen. Based on the flux-calibrated (but low- 
resolution) B&C spectrum, we can derive the line luminosities 
of the [Oi] lines. Zickgraf et al. (119861) determined an interstellar 
extinction towards S 65 of E(B - V) = 0.18 ± 0.02 (see Table[IJ. 
We used this value and applied the interstellar extinction curve of 
the SMC of Prevot et al. (1984) to deredden the B&C spectrum. 
With the calibrated dereddened continuum and a distance mod- 
ulus of /i = 18.93 + 0.02 for the SMC (Kefler & Wood l2006] l. 
we then calculated the observed luminosities in the [Oi] lines 
(see TableO from the measured equivalent widths in the FEROS 
spectrum. The errors in line luminosities are about 10%, follow- 
ing from the uncertainties in equivalent width measurements, in- 
terstellar extinction, and distance. 



[A] [kms-'] [kms-'] [ergs"'] 



5577 
6300 
6364 



44 ±2 
24 ±2 
24 ±2 



56: 

43: 
42: 



8.28 X 10"" 
8.46 X 10"* 
2.93 X 10"* 



The major advantage in using forbidden emission lines is that 
because they are optically thin, their line profiles contain the full 
velocity information of their formation region. The observed line 
profiles of the [Oi] A 6300 and 6364 lines are more or less iden- 
tical. This is clear, because both lines arise from the same upper 
level, which means that they are formed in the same region in 
terms of density and temperature. While both lines show only 
a slight indication of a double-peaked profile, this double-peak 
structure is more pronounced in the [Oi] /1 5577 line, and the 
measured peak velocities and FWHM values of all three lines 
are listed in Table[3] This latter line arises from a higher energy 
level. Because the levels from which the forbidden lines origi- 
nate are purely collisionally excited, the [Oi] A 5577 line must 
be formed in regions of a higher density and/or higher tempera- 
ture. This could indicate that the /1 5577 line originates from re- 
gions closer to the star. In combination with the broader FWHM 
as well as the broader peak seperation, which mirror a higher 
velocity within the A 5577 line-forming region, this qualitative 
analysis might already hint at a Keplerian disk rather than a pure 
outflowing disk. 

However, based purely on the appearance of a double-peaked 
structure of the line profiles, it is not possible to discriminate 
whether the [Oi] line-forming region is located in a Keplerian 
rotating disk or in an outflowing disk, because both scenarios, 
when viewed edge-on, result in identical, double-peaked line 
profiles. Reliable conclusions about the real kinematical nature 
of the [Oi] line-forming region can thus only be drawn from a 
proper and simultaneous modeling of both line profiles, from 
which the kinematical information is taken and the line lumi- 
nosities, delivering the information about the density and tem- 
perature. This is done below for the two competing scenarios: 
the outflowing and the Keplerian rotating disk model. 

3.3.1. Outflowing disk model 

Not much is known in the literature concerning the structure and 
shape of the disks around B[e] supergiants, and there is an on- 
going debate whether they are Keplerian rotating, or outflowing 
in nature (see, e.g.. Porter l2003l Kraus et al. 120071) . From a sta- 
tistical point of view (based on low number statistics though), it 
seems that the disks must be rather thick, with an opening an- 
gle of about a ^ 20° (Zickgra f 1989) , which is 4-5 times larger 
than the opening angle of classical Be stars (see, e.g.. Porter & 
Rivinius 2003]) . 

Another key question is related to the ionization structure 
in the disk. Oxygen has about the same ionization potential as 
hydrogen. The detection of emission lines from neutral oxygen 
thus means that these lines must be generated within a region in 
which hydrogen is predominantly neutral as well. Free electrons 
are (besides the less efficient, but nevertheless important, neutral 
hydrogen atoms) the main collision partners to excite the lowest 
energy levels in neutral oxygen, from which the forbidden emis- 
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sion lines arise. We follow the approach of Kraus et al. (I2007I I 
and introduce the parameter ^e as the ionization fraction, de- 
fined by «e - qetiii, with ^e = ^Metals + ?H+, and the electron and 
hydrogen particle densities, n^ and n-a, respectively. Assuming 
that all elements with an ionization potential lower than that of 
hydrogen and oxygen are fully ionized delivers an upper limit of 
^Metals < 9.5 X 10"^ For this, we used an SMC metallicity of 20% 
solar (see, e.g., Dufton et al. 2005), and solar abundance values 
for the metals from Grevesse & Sauval (119981 ). This very low 
value obtained from the metals means that if hydrogen is ion- 
ized by only 1%, it still provides the dominant quantity of free 
electrons. Because we do not know the amounts of ionized hy- 
drogen and ionized metals within the [Oi] line-forming region, 
we keep g'e as a free parameter at first. 

An additional important parameter is the electron tempera- 
ture in the [Oi] line-forming region. It should be lower than the 
hydrogen ionization temperature, and higher than the molecule 
dissociation threshold, below which atomic oxygen is bound into 
CO molecules. While the former restriction defines the maxi- 
mum temperature at about lOOOOK, the latter corresponds to a 
minimum temperature of roughly 5000 K, which is the dissocia- 
tion temperature of CO. 

To calculate the population of the five lowest energy levels 
in atomic oxygen, we need to take into account collisions from 
both free electrons and neutral hydrogen atoms. The collision 
parameters for excitation by neutral hydrogen atoms are taken 
from Storzer & Hollenbach (2000), those for excitation by free 
electrons from Mendoza ( 1983 ), and atomic parameters are from 
Wiese et al. (1966) and Kafatos & Lynch ( 119801 ). We calculate 
the level population by solving the statistical equilibrium equa- 
tions in a 5-level atom. Forbidden emission lines are optically 
thin so that no complicated radiation transfer needs to be solved. 
Instead, the emissivity at any location r in the disk is simply 
given by j„,,„(r) = hv„j„n„{r)A„^„, where n„{r) is the level popu- 
lation of the upper level, and A„ „, is the Einstein coefficient of 
spontaneous emission. The line luminosities finally follow from 
integration of the emissivities over the emitting disk volume. 

Next, we can restrict the valid ranges in electron tempera- 
ture and ionization fraction, q^, within the [Oi] line-forming re- 
gions in the disk of S 65. Because two of the observed [Oi] lines 
(/16300 and /16364) originate from the same upper level, their 
line ratio is determined by pure quantum mechanics and does 
not depend on density and temperature. The A5577 line, how- 
ever, originates from a higher energy level. A sufficiently strong 
population of this level thus severely depends on the density of 
the collision partners as well as on temperature. Any ratio with 
this line is therefore an ideal tracer for the density and tempera- 
ture, and we concentrate on the /16300//15577 line ratio below to 
constrain these two parameters. 

In the top panel of Fig.]?] we show the /16300//15577 line 
ratio as a function of the ionization fraction, q^, and for dif- 
ferent (constant) values of the electron temperature. The calcu- 
lations performed with the outflowing disk model request the 
knowledge of the mass flux. For simplicity, we assume a lat- 
itude independent mass flux within the disk-forming wind re- 
gion, i.e., F^{0) - Fn,,Disk - const, for all latitudes in the range 
[90° - (a/2)] <0 < [90° + (a/2)]. We also assume the wind to 
have a constant ouflow velocity, iiout, delivering the density pa- 
rameter, /^m.Disk/fout, as One free input parameter. With this den- 
sity parameter we can calculate the hydrogen density at the inner 
edge of the disk (here: at the stellar surface) and consequently 
also the radial density distribution nu(r) ~ (^m.Disk/i'out) f^^ re- 
sulting from the application of the equation of mass continuity. 
Due to the radially decreasing density, the [Oi] line luminosi- 
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Fig. 7. Top panel: Theoretical line ratios as function of q^ at 
fixed temperature. The intersection with the observed value (dot- 
ted horizontal line) defines the valid (^e, T^) range (arrows). 
Bottom panel: Density parameter as function of q^ at fixed tem- 
perature. The valid (^e, ^e) range (dotted vertical lines) delivers 
the range in density parameters (arrows). 



ties saturate at a certain distance. This saturation value should 
agree with the observed line luminosities. We thus varied the 
density parameter, /^m,Disk/fout, for a fixed pair of (T^, q^) val- 
ues, until the observed /16300 line luminosity was reproduced. 
The line ratio with the A5511 line resulting from these calcula- 
tions defines one single point in the diagram shown in the top 
panel of Fig.|7] This procedure was repeated until a sufficiently 
large coverage of the (T^, q^) parameter space was obtained. We 
then compared the theoretical results to the observed line ra- 
tio, which is shown as the dotted line. Based on the discussion 
above concerning the valid range in temperatures, the curve for 
Tf. = 5000 K can be considered as the lower temperature limit, 
resulting in a strict upper limit for the ionization fraction in the 
disk of S65 of g'e.max = 4.2 X 10"^. This means that the [Oi] 
line-forming region must indeed be neutral in hydrogen, with 
q\i+ ~ ?e,max = 4.2 X 10"^. For a temperature of 7000 K, the 
observed line ratio cannot be reproduced. The maximum tem- 
perature in the [Oi] line-forming region must thus lie between 
6500 K and 7000 K. Since the plot in the top panel of Fig.|2]in- 
dicates that a higher temperature would imply a lower degree 
of ionization, it is reasonable to assume that the most plausible 
range in temperature is Tg ^ 5500 . . . 6500 K, delivering a disk 
ionization fraction of g'e - (1 . . .4) x 10"^. 

In the bottom panel of Fig. [T] we plot the density parameter, 
used for the line luminosity fit of the /16300 line and the resulting 
line ratio as shown in the top panel, as a function of q^ and for 
the different temperature values. In the part of the diagram where 
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the ionization fraction is high, the density parameter shows (for a 
given temperature) a power-law dependence with q^. This means 
that for decreasing ionization fraction the input density has to be 
increased to compensate for the otherwise decreasing number of 
free electrons. On the other hand, below an ionization fraction 
of a few times 10"^, the input density parameters stagnate. This 
means that for ionization fractions q^ % lO""* collisions with 
neutral H atoms start to become important. For further decreas- 
ing ionization fraction they even dominate the level population 
process, because at very low ionization the neutral particles are 
finally numerous enough for efficient collisional level popula- 
tion. From the narrow range of valid temperatures and ioniza- 
tion fraction values, the range in density parameter is confined 
to F™,Disk/i^out ^ (1 ... 2) X 10-"' g cm-l 

These calculations, even though they are performed under 
several simplifying assumptions, show that the [Oi] line-forming 
region is well constrained in temperature and density. And the 
increase of the individual line luminosities with distance from 
the star and their saturation is shown for all three [Oi] lines 
in the upper panel of Fig. [8] These computations exemplify the 
case of a constant temperature of 6000 K and resulting values of 
2 .49 X 1 0"^ and 1.16xl0-"^g cm"-' for the ionization fraction and 
density parameter, respectively. The arrows to the right indicate 
the observed line luminosity values with their errors indicated 
by the vertical bars. 

The luminosity saturation of the A5511 line occurs at a dis- 
tance of about 20 - 30 /?t, and therefore rather close to the star, 
while saturation in the luminosities of the other two lines hap- 
pens only at distances beyond ~ 400/?,, i.e., at least 10 times 
further out and therefore in regions in which the density is at 
least 100 times lower To check the influence of the disk temper- 
ature on these results we also calculated the line luminosities for 
a disk with 6500 K and 5500 K. A lower (higher) temperature 
requires slightly higher (lower) values for the ionization fraction 
and density parameter (see Fig.|7]l. The corresponding curves for 
the increases in line luminosities are, however, hardly affected, 
and the saturations of the [Oi] lines happen at the same distances. 
What remains to check is whether this outflowing disk scenario 
is a reasonable model to explain the observed line profiles. 

Because the lines are formed in different disk regions, their 
profiles mirror the kinematics at these different distances from 
the star. The double-peaked profile of the A5511 line might be 
interpreted with an outflow velocity of about Wout - 22kms"', 
obtained from the separation between the blue edge of the blue 
and the red edge of the red peak (TableO. The broader wings 
of the line require an additional Gaussian-shaped velocity com- 
ponent, i.e., additional to the thermal velocity plus spectral res- 
olution, which result in a Gaussian velocity component of about 
6 km s"' . This additional Gaussian velocity could be ascribed to 
some turbulent motion of the gas and is found to be on the order 
of ytmb - 8 kms"'. The resulting fit to the A5511 line is shown 
as the dashed line in the bottom panel of Fig. [8] 

For the other two [Oi] fines (/16300 and /16364) this model 
does not hold (dotted lines in the lower panel of Fig.[8]l. Instead, 
the observed profiles are much narrower and the peak-separation 
is less pronounced than in the /15577 line. Nevertheless, we find 
good fits to these line profiles, too, if we assume an outflow ve- 
locity of only 16km s', but a turbulent velocity of llkms ' 
(dashed lines in the lower panel of Fig. [8). This would imply 
that because these lines are generated much further out, the out- 
flow velocity of the disk would have slightly slown down, with 
a simultaneous slight increase of the turbulent velocity. 
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Fig. 8. Top panel: Predicted [Oi] line luminosities from within a 
radius r in the outflowing disk scenario. The arrows to the right 
indicate the observed values. Bottom panel: Fits to the line pro- 
files. The dotted profiles fitted to the /16300 and /16364 lines are 
obtained for the same constant outflow velocity as the A5511 
line; the dashed lines in all profiles show our best fit. 



3.3.2. Keplerian rotating disk or ring 

The decrease in peak separation seen in the [Oi] lines from a 
higher value at higher densities, i.e., closer to the star, to a lower 
value further out, could be a clear indication for Keplerian rota- 
tion. Therefore we now calculate the [Oi] line luminosities and 
the profiles with a Keplerian rotating disk model. To translate the 
rotation velocities as inferred from the peak separations of the 
line profiles into radial distances, we need to know the current 
mass of the star Its effective temperature and luminosity, when 
compared to evolutionary tracks of rotating (Meynet & Maeder 
2005 ) and non-rotating stars at SMC metallicity (Charbonnel et 
al. .l993j , place S 65 within the HR diagram to a region in agree- 
ment with a post-main sequence evolutionary phase of a star with 
an initial mass in the range 35^0 M©, and we assume a current 
mass of ~ 35 M©. 

The peak separations of the /16300 and /16364 lines de- 
liver a minimum rotation velocity of about 5 ... 6 km s"' , which 
sets the outer edge of the [Oi] line emission region to about 
Rom - 3000 R*. The A5511 line, which is created much closer 
to the star, thus delivers the velocity information from the inner 
parts of the Keplerian rotating disk. Due to the extended region 
over which the [Oi] line is generated, its peak separation does not 
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Fig. 9. Top panel: Same as Fig.[8]but for the Keplerian rotating 
disk scenario. Bottom panel: Fits to the line profiles. The dotted 
profile fitted to the A5511 line is for a purely rotating disk model; 
the dashed lines in all profiles account for our best-fit model. 



reflect the rotation velocity at the inner edge of the disk, but the 
velocity at the distance where the bulk of emission in this line is 
created. The distance at which a rotation velocity of 22 km s"' is 
achieved is at about 170R*. The FW//M of about 56 km s"' im- 
plies a rotation velocity of 28 km s"' and a distance of 105 R* at 
which already some substantial fraction of the emission must be 
generated. The real inner edge of the disk must be slightly closer 
to the star though, and our test calculations revealed a most likely 
distance of Rm - 90 R* delivering a maximum rotation velocity 
to the profile of the /I5577 line of Urot.max - 30km s'. 

While the velocity distribution in a Keplerian rotating disk is 
well known, the situation is far less clear for the radial density 
distribution within this disk. For simplicity we adopt a r^^ be- 
havior of the density distribution. A justification for this choice 
might be given by the request that the disk must have been 
formed from an outflow or ejection process obeying the mass 
continuity relation. In addition, a r^^ density distribution is also 
known to settle over large regions within the so-called wind- 
wind interaction regions (e.g., Steffen & Schonberner I2003I) . 
However, different power-law density distributions, as, e.g., 
found in Keplerian disks around pre-main sequence stars, can- 
not be excluded. 

For a first guess of the ionization fraction and density pa- 
rameter at the inner edge of the disk, we calculate the line lu- 
minosities and the /16300//15577 line ratio for the same range 
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Fig. 10. Expected strength and shape of the [On] /i7319 line in 
the Keplerian rotating disk scenario. 



of constant temperatures (5000 K to 7000 K) and a disk extend- 
ing from ~ 90 R* to ~ 3000 R,. Proceeding as in Sect. l3.3.ri we 
find the following plausible ranges for the disk parameters: q^ ^ 
(1 .3 ... 3.6) X 10-3 and F„,Disk/fout ^ (0.5 . . . 1 .0) x 10"'° g cm-^ 
or, when expressed in terms of the density, p, at the inner edge 
of the disk,p(/?in) ^ (0.6 ... 1.2) x lO^'^^gcm-^ 

The results of the line luminosity calculations performed for 
a constant temperature of Te - 6000 K are shown in the top 
panel of Fig.|9] They indicate that for a Keplerian rotating disk 
the ionization fraction, q^^, must be higher by about a factor of 
100 compared to the one found for the outflowing disk model. 
This trend is clear, because the line formation is now forced to 
happen at much larger distances from the star, where the Oi den- 
sity is generally much lower Therefore, a higher number of col- 
lision partners (i.e., free electrons) is necessary to produce the 
same line luminosity. Nevertheless, even with this much higher 
ionization fraction, the disk can still be considered as being neu- 
tral in hydrogen (^h+ < 2 x lO""*). 

The calculated profile fits to the [Oi] lines, resulting from 
our Keplerian disk model, are shown in the bottom panel of 
Fig.|9] It is obvious that the /16300 and /16364 lines can be fitted 
quite well, as shown by our dashed-line fits to these two lines. 
However, the same model fails (indicated by the dotted line) to 
reproduce the A5511 line profile. Although the Keplerian disk 
model results in a slightly double-peaked profile for that line as 
well, its shape is nevertheless significantly diff'erent from the ob- 
served one. Especially the inter-peak strength is strongly over- 
estimated. To obtain a better reproduction of the A5577 line pro- 
file, we thus have to assume that the material at the inner edge 
of the disk is not purely rotating, but an additional kinematical 
component, resulting in a double-peaked line profile, must exist 
as well. Because our investigations in the previous section have 
shown that an equatorial outflow produces double-peaked line 
profiles, we can speculate that the disk material is not revolv- 
ing on stable Keplerian but rather on slightly outwards spiraling 
orbits. 

Testing this hypothesis revealed that to reproduce the /15577 
line profile we need to include an additional outflow veloc- 
ity at the inner edge of the [Oi] line forming region of Uom = 
(9.0 + 0.5)kms '. The resulting line profile is shown as the 
dashed line fit to the /15577 line in the lower panel of Fig.|9] 
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Table 4. Best-fit parameters of the two competing scenarios for 
S65. 



Table 5. Velocities obtained from the peak separation and the 
FWHM values, and equivalent widths of the [Can] emission 
lines. 



Parameter 


Outflowing disk 


Keplerian disk 


A 


"peaks 


"FWHM 


EW 




'^in "out 


^in 


'^oul 


[A] 
7291 


[kms-'] 
42.0 ± 2 


[kms-'] 

55.2 ± 2 


[A] 


r[R,] 


1 >400 


-90 


-3000 


3.04 ±0.21 


youtflow [kms '] 


22 ±1.0 16 ±1.0 


9 ±0.5 


^2.0 


7324 


38.6 ±2 


56.8 ± 2 


2.09 ±0.11 


Wrot [kms"'] 


— — 


30 ± 1.0 


5 ±0.5 










o,„,b [kms-i] 


8 ±0.5 11 ±0.5 


— 


— 










n^(r) [cm-3] 


1.2 X 10' < 7.6 X 103 


8.2 X 10* 


7.4 X 10^ 











On the other hand, the profiles of the lines AA6300, 6364, which 
are formed at a much larger distance, are not compatible with 
such an outflow component. In principle no outflow is needed 
for these lines, but an outflow component of Uout ^ 2 kms-' 
would still reproduce the profiles reasonably well, indicating a 
slow-down of the outflow. 

If the observed [Oi] lines result indeed from a Keplerian (or 
quasi-Keplerian) rotating disk, the question arises whether this 
disk extends down to the stellar surface, or whether the mate- 
rial is detached from the stellar surface. If the disk extends down 
to (or close to) the stellar surface, the disk parts within 90/?», 
from which no contribution to the [Oi] emission is seen, should 
consequently be ionized. Because ionized oxygen is also known 
to have forbidden emission lines in the optical spectral range, 
we calculated the expected line luminosity and profile for the 
strongest of these lines, which is usually the [On] /17319 line, 
arising in the inner disk parts that extends from the stellar sur- 
face out to 90 Rf The resulting line profile and intensity are then 
compared to the weakest of our [Oi] lines, i.e., the ASSll Une. 
This comparison is shown in Fig. [10] The [On] line turns out to 
be much broader (due to the higher rotation velocities closer to 
the star) and about five times more intense than the /15577 line. 
Therefore it should be clearly detectable. But our spectra do not 
show indications for the presence of any [On] line. This could 
indeed indicate that if the Keplerian rotating disk scenario is the 
correct one, the material must be detached from the star. 



4. Discussion 

4.1. Outflowing versus Keplerian rotating disk 

So far, modeling of the [Oi] line luminosities and profiles alone 
did not allow us to unambiguously distinguish between the out- 
flowing and the Keplerian rotating disk model. Both can provide 
a plausible scenario for the neutral material around S 65, and we 
summarized the best-fit model parameters for each in Table|4] 
The question thus arises whether we can find additional, pre- 
dominantely kinematical hints in our data that would allow us to 
favor one of the disk models. 

Besides the [Oi] lines, our spectra display clearly double- 
peaked line profiles of the [Can] lines (see Fig.|4|. The measured 
velocities and equivalent widths are listed in Table|5] The ob- 
served line ratio [Can] /17291/.17324 of 1.45 + 0.13 does not 
give a conclusive hint for the formation region of these lines, 
but agrees more with the formation in the low-density limit 
(1.495) rather than in the high-density limit (1.535). These val- 
ues have been determined and discussed in detail by Hartigan et 
al. (12004 ). and the critical density is ~ 5 x 10^ cm--*. 

The measured peak and FWHM velocities are very similar 
to the values found for the [Oi] A5511 line, indicating that the 
[Can] lines originate from the same region. In both disk scenar- 



ios the electron density in the [Oi] line formation regions are the 
same, i.e., in the region of the [Oi] A5577 line formation it ranges 
from He - 8.0 x 10^ ... 1.5 x 10^ cm--', while it is much lower 
within the [Oi] AA63QQ, 6364 line formation region, extending 
over Me - 5.0 x 10^ . . . 7.5 x lO"' cnT^. Consequently, the [Can] 
lines are indeed formed below the critical density as expected 
from the measured line ratio. Although the kinematics derived 
from the [Can] lines cannot help to distinguish between the out- 
flowing and the Keplerian rotating disk we emphasize that the in- 
tense and clearly double-peaked [Can] lines seen in the spectrum 
of S 65 provide an additional strong indication for the existence 
of a gaseous disk. 

Nevertheless, we think that the Keplerian rotating disk is the 
more reasonable scenario. Supporting kinematical arguments in 
favor of the Keplerian disk are given by: (i) the symmetric, un- 
shifted forbidden emission lines and (ii) the sharp absorption 
components. 

We measured the FWHM velocities of all unblended for- 
bidden emission lines. The velocities of all lines are very sim- 
ilar, resulting in a mean value of (37.09 + 4.16)kms-'. This 
value is slightly lower than the FWHM velocities of the [Oi] 
AA6300, 6364 lines, indicating that the bulk of forbidden emis- 
sion lines originates in regions near but beyond (i.e. at larger 
distances and hence lower velocities than) the [Oi] lines. This 
scenario is difficult to explain with an outflowing disk model, 
especially because that model requires an increase in turbulent 
velocity with distance from the star. On the other hand, it is logi- 
cal to explain the profiles of the forbidden emission lines with the 
Keplerian disk model. The symmetric line profile with no peak 
separation means that the emission is formed over a large disk, 
as for the [Oi] AA63QQ, 6364 lines, but extending to distances at 
which the rotation is too small (< 5 kms-') for the peaks to be 
still separated. 

Even stronger support for the Keplerian disk scenario is 
given by the numerous shell lines, i.e., the emission lines with 
a sharp absorption component reaching below the stellar undis- 
turbed flux. As for classical Be shell stars (e.g.. Porter & Rivinius 
2003), we can use these sharp absorption components to ex- 
tract the information on the radial velocity structure, because 
these lines involve absorption of the photospheric flux through- 
out the whole disk. The wavelength displacements of the absorp- 
tion components with respect to the line center of the emission 
line together with the measured widths of the absorption compo- 
nents deliver information about the kinematics of the absorbing 
medium, i.e., the kinematics of the disk. 

For all the unblended emission lines with pronounced ab- 
sorption components we measured the central wavelengths and 
the FWHM values of their absorption components. The cen- 
tral wavelengths of the absorption components were corrected 
for the systemic velocity. The velocities resulting from their 
FWHM values are plotted in Fig. [TT] versus their velocity offsets 
of their central wavelengths, defined by the difference Uobs - fsys- 
From this plot we can draw the following conlusions: Firstly, 
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all absorption components are blueshifted, meaning that those 
disk parts in which the absorption happens show a slow outflow 
ranging from to about 5kms"'. These low outflow velocities 
agree well with those found for the Keplerian disk scenario (see 
Table|4]l. Secondly, the FWHM velocities of the absorption com- 
ponents range from about 6 to about 13 km s"' . Interpreting them 
with the rotation velocity of the disk, the bulk of the absorption 
happens at rotation velocities between about 3 and 6.5 km s"'. 
In addition, there seems to be a slight trend (but this trend might 
also be spurious and caused by a few outliers) towards slower 
outflow for lower rotation velocities as indicated by the dashed 
line, which represents the linear regression calculated for the 
data points. In that respect it is also interesting to note that this 
trend is also seen when comparing the results for Feii with those 
of the (though rare) Fei lines. The latter are concentrated in the 
right part of the plot, i.e., at the region with lower outflow and 
lower rotation velocities, which is logical, because we expect an 
outwards decreasing ionization structure in the disk. 

A further point that needs to be discussed is the inclination 
of the system. All our conclusions presented so far are based on 
an assumed edge-on orientation of the system. As discussed in 
Sect. l3.2l the derived projected rotation velocity is not conclusive 
and can only be considered as a lower limit to the real projected 
velocity. Even if we request the star to rotate at its critical limit 
(see Sect. l4. 21 below), we cannot claim that the inclination then 
has to follow from / = arcsin((;i-ot/i'crit) - arcsin(0.75) = 48.6°, 
but it can have any value between 48.6° and 90°. On the other 
hand, as mentioned in Sect. 13. II the numerous shell lines speak 
in favor of a system viewed more or less edge-on. With a to- 
tal opening angle for B[e] supergiant stars' disks of ~ 20° (see 
Sect. 13. 3. Ti l, a minimum inclination of about ~ 80° might still 
be possible for the line-of-sight to pass through the disk mate- 
rial. In this case, the difference in velocity compared to a com- 
pletely edge-on orientation is negligible. Yet we have to admit 
that this argument only holds if the disk maintains its constant 
opening angle, i.e. if its outer parts are not flared. So far noth- 
ing is known about the structure and shape (and especially about 
a possible flaring) of the disks around B[e] supergiants at large 
distances. Therefore we consider the close to edge-on orientation 
(i.e. / <: 80°) used in our model as a plausible scenario. 

Although the above summarized arguments can only be re- 
garded as qualitative indications rather than real proofs, they 
nevertheless seem to speak in favor of a Keplerian disk scenario. 
What remains to be discussed is therefore the problem of the 
origin of such a disk. 

4.2. Possible disk formation mechanism 

The source S 65 is in an evolved evolutionary phase of an 
early-type progenitor star, with an initial mass of 35 . . .4OM0. 
Hydrodynamic calculations for a 35 M© star by Garcia-Segura 
et al. (11996b) have shown that due to its high mass-loss rate and 
wind speed during its main-sequence evolution, the stellar wind 
deposits sufficient energy into the surrounding medium to sweap 
up the circumstellar matter into a so-called interstellar (or main- 
sequence) bubble. When entering the supergiant phase, this bub- 
ble has grown in size up to a radius of about 35 ... 40 pc (Garcia- 
Segura et al. 1996b). Consequently, the disks seen around B[e] 
supergiants with sizes of only a few hundred to a few thousand 
AU cannot be pre-main sequence in origin. 

In single-star evolution, the disk-formation mechanism in 
evolved massive stars is not well understood (for a review see, 
e.g., Owocki 2006). However, as first noted by Struve ( |1931| l 
rapid (critical) stellar rotation seems to play the key role. Studies 
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Fig. 11. FWHM velocities versus central peak velocities of the 
sharp absorption components. Errors in velocities are indicated 
by the cross in the lower right corner. The dashed line indicates 
the result from a Hnear regression. 



on disk formation around the most rapidly rotating stars known, 
i.e., the classical Be stars, revealed that their disks are most 
probably driven by viscosity (see, e.g.. Porter 1999, Okazaki 
2001. Jones et al. .2008J . resulting in quasi-Keplerian rotating 
disks with a slow (few km s"') equatorial outflow component, in 
agreement with observations (Clark et al. 2003). The star S 65 
also seems to be rotating (close to) critical. Consequently, only a 
slight perturbance, perhaps in form of stellar pulsations (Owocki 
i2006i ), would be sufficient to launch matter into orbit and to ini- 
tiate a viscously driven disk. The kinematics found for the disk 
of S 65 seem to agree with a quasi-Keplerian disk with a low 
outflow velocity. 

Other possibilities for disk formation would usually require 
a close companion, which cannot efficiently accrete the mate- 
rial transferred from the primary so that an accretion disk would 
form. The best known example for this mass transfer and disk- 
formation process in a massive evolved system is certainly 77 Car 
(e.g. Soker l2003]) . But though studied for decades, none of the 
B[e] supergiants has so far been reported to have a close com- 
panion (e.g. Kraus et al. 1201 Oi l. In addition, neither the spectral 
energy distribution of S 65 nor our high-quality spectra, which 
have been taken on dififerent dates covering observing periods 
as long as 1 1 years and as short as 24 hours, show any indica- 
tion for a companion. Also, our spectra show no evidence for 
ongoing accretion. Although the existence of a companion can- 
not definitely be ruled out, we consider it as the less probable 
scenario for the disk formation around S 65. 

Whether the disk around S 65 is indeed driven viscously, 
as suggested by the observed (quasi-)Keplerian rotation and the 
slow outflow component, needs to be studied in a future work. 
Also, whether the disk indeed slows down from inside out as 
suggested by our data, is an interesting question that needs to be 
investigated in much more detail than can be done here. 
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4.3. The connection with LBVs 

Another still unsolved problem is related to the question whether 
B[e] supergiants are a special evolutionary phase in single-star 
evolution through which only massive stars with very specific 
initial conditions will evolve, or whether all massive stars have 
to go through the B[e] supergiant phase. And if the latter is true, 
then what are the direct progenitors/descendants of the B[e] su- 
pergiant phase ? 

B[e] supergiants are often discussed in connection with 
Luminous Blue Variables (LBVs). This is certainly because 
LBVs also show an indication for a disk or ring structure of 
circumstellar material. Although a bipolar or multipolar struc- 
ture seen in several LBVs like in 77 Car can be caused by binary 
interaction (Soker 2003 ), the majority of LBVs have not been re- 
ported to have a close companion. Instead, their often seen ring- 
shaped nebulae (e.g., Nota et al. 119951) can easily be explained 
by wind- wind interaction in single-star evolution (Garcia-Segura 
et al. ll996ab . Morphological studies of LBV nebulae imply that 
the slow and dense wind arises mainly from the equatorial re- 
gions (Nota et al. 1995) . And as for the B[e] supergiants, rapid 
stellar rotation is thought to play a major role in triggering this 
non-spherical mass-loss. 

The second common characteristic is related to the outflow 
velocity. Weis (120031) studied the kinematics of LBV nebulae 
and found that their expansion velocities seem to decrease with 
metallicity. The LBV nebulae in the SMC are thus expected to 
have expansion velocities of less than 10 km s"' ; hence S 65 with 
its outflow component of a few kms ' would fit in perfectly. 
With regard to this, we add that for the B[e] supergiant R 126 
in the LMC a disk outflow velocity of ~ 1 1 .5 km s"' was deter- 
mined (Kraus et al. 120071 ). in good agreement with the outflow 
velocities of LBV nebulae in the LMC. We also find that the gas 
masses in the disks around S65 (^ 1.5 x 10"^ Mq) and R126 
(^ 6 X 10"^ Mq) fit nicely to the range of gas masses found for 
LBV nebulae (Nota et al. lT995l Weis l2003T ). 

The last common characteristic refers to the stellar rotation 
speed. Groh et al. (2009 ) recently found for two of the Galactic 
LBVs, AG Car and HR Car, rotation velocities during their vi- 
sual minimum phase in excess of 85% of their critical velocity. 
Hence Groh et al. (120091 ) defined an LBV minimum-instability 
strip, which indicates the position in the HR diagram at which 
critical rotation for LBVs is expected. This is shown in Fig. [12] 
To the left of the instability strip is the physically unstable re- 
gion for LBVs, labeled with Urot/f'crit > 1- We also included the 
position of S 65 into this plot. Its location is in the valid range 
for LBVs, but also very close to the instability strip. 

In summary, the rapid rotation, the close vicinity in the HR 
diagram to the LBVs and the LBV minimum instability strip, 
as well as the similar disk mass and outflow velocity derived 
for S 65 might all hint towards a connection between LBVs and 
S 65 (or perhaps B[e] supergiants in general). The most logical 
conclusion is thus that S 65 is currently in the LBV visible mini- 
mum phase, after experiencing one (or more) eruptions, in which 
a large amount of material was ejected from the rapidly rotat- 
ing star, predominantly in equatorial direction, which formed the 
Keplerian ring or disk-like structure. If this scenario is correct, 
then S 65 could be in a state just before the new, fast wind sets 
in, which will compress the disk/ring material into the typical 
LBV nebula. The signatures of Keplerian rotation, as still seen 
in the disk of S 65, might disappear during this wind compres- 
sion process, because the material will be radially accelerated 
outwards. 
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Fig. 12. Location of S 65 within the HR diagram and with re- 
spect to two rapidly rotating Galactic LBVs close to the LBV 
minimum instability strip (after Groh et al. 120091 ). For details 

see text. 



We cannot conclusively answer whether S 65 indeed belongs 
to the class of LBVs, especially because a star is called an LBV 
only after a giant eruption and the typical S Dor variability has 
been recorded. Nevertheless, there is ample evidence that S 65 
shares several of the major LBV characteristics. 



5. Conclusions 

We studied the circumstellar material of the rapidly rotating 
SMC B[e] supergiant star S65. The sharp absorption compo- 
nents seen in its numerous emission lines indicate a high-density 
disk viewed almost edge-on. Surprisingly, the spectrum also dis- 
plays two very intense and clearly double-peaked [Can] lines, 
supporting the high-density gaseous disk hypothesis. Based on 
the modeling of both, the observed line luminosities of the [Oi] 
lines and their line profiles, we found that the disk around the 
B[e] supergiant S 65 must be predominantly neutral in hydrogen. 
The total amount of ionized hydrogen is found to be less than 
~ 0.1%, and the total amount of gas traced within the [Oi] line- 
forming region is ~ 1 .5 x 10"^ M©. In combination with the kine- 
matical information provided by the numerous lines, our model 
results favor the interpretation of a detached (quasi-)Keplerian 
rotating disk rather than an outflowing disk. The location of S 65 
in the HR diagram places the star close to the recently defined 
LBV minimum instability strip. Together with its high rotation 
velocity close to the critical limit and its large amount of slowly 
expanding circumstellar material similar to that seen for LBV 
nebulae, it seems that S 65 might be the link between B[e] su- 
pergiants and LBVs. 
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